Time series technical analysis via new fast 
estimation methods: a preliminary study in 
mathematical finance 
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Abstract: New fast estimation methods stemming from control theory lead to a fresh look at 
time series, which bears some resemblance to "technical analysis". The results are applied to 
a typical object of financial engineering, namely the forecast of foreign exchange rates, via a 
"model-free" setting, i.e., via repeated identifications of low order linear difference equations on 
sHding short time windows. Several convincing computer simulations, including the prediction 
of the position and of the volatility with respect to the forecasted trendline, are provided. Z- 
transform and differential algebra are the main mathematical tools. 
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1. INTRODUCTION 



1.1 Motivations 



Recent advances in estimation and identification (see, 
e.g., (Fliess & Sira-Ramirez, 2003; Fliess & Sira-Ramirez, 
2004; FHess, Join & Sira-Ramirez, 2004; FHess & Sira- 
Ramirez, 2008; FHess, Join & Sira-Ramirez, 2008) and the 
references therein) stemming from mathematical control 
theory may be summarized by the two following facts: 

• Their algebraic nature permits to derive exact non- 
asymptotic formulae for obtaining the unknown quan- 
tities in real time. 

• There is no need to know the statistical properties of 
the corrupting noises. 

Those techniques have already been applied in many 
concrete situations, including signal processing (see the 
references in (Fliess, 2008)). Their recent and successful 
extension to discrete-time Hnear control systems (FHess, 
Fuchshumer, Schoberl, Schlacher & Sira-Ramirez, 2008) 
has prompted us to study their relevance to financial time 
series. 

Remark 1.1. The relationship between time series analysis 
and control theory is well documented (see, e.g., (Box, 
Jenkins & Reinsel, 1994; Gourieroux & Monfort, 1995; 
Hamilton, 1994) and the references therein). Our view- 
point seems nevertheless to be quite new when compared 
to the existing literature. 

Remark 1.2. The title of this communication is due to 
its obvious connection with some aspects of technical 



analysis, or charting (see, e.g., (Aronson, 2007; Bechu, 
Bertrand & Nebenzahl, 2008; Kaufman, 2005; Kirkpatrick 
& Dahlquist, 2006; Murphy, 1999) and the references 
therein) , which is widely used among traders and financial 
professionals. ^ 

1.2 Linear difference equations 

Consider the univariate time series {x{t) | t G N}: x{t) 
is not regarded here as a stochastic process like in the 
familiar ARMA and ARIMA models but is supposed to 
satisfy "approximatively" a linear difference equation 

x{t + n) — aix{t -\- n — 1) — ■ ■ ■ — anx{t) = (1) 

where ai,...,a„ G M. Introduce as in digital signal 
processing the additive decomposition 

X{t) = a;trendline(i) + '^(i) (2) 

where 

• a;trendiine(0 IS the trcndUnc^ which satisfies Eq. (1) 
exactly; 

• the additive "noise" i^{t) is the mismatch between the 
real data and the trendline. 



Thus 

x{t -\- n) ~ aix{t + n — 1) — • • • - 
where 

e(t) = v{t + n) — aivit + ?i — 1) 



a„x(i) = e{t) (3) 



anv{t) (4) 



Technical analysis is often severely criticized in the academic 
world and among the practitioners of mathematical finance (see, e.g., 
(Paulos, 2003)). 

^ Compare, e.g., with (Durlauf & Philips, 



We only assume that the "ergodic mean" of v{t) is 0, i.e 
v{0) + + ■ ■ ■ + iy{N) 



hm L_21_^ = 

AT^+oo N +1 

It means that, V i G N, the moving average 



MA,,jv(i) = 



i^{t)+iy{t + l) + --- + v{t + N) 



(5) 



(6) 



N + 1 

is close to if iV is large enough. It follows from Eq. (4) 
that e{t) also satisfies the properties (5) and (6). Most of 
the stochastic processes, like finite linear combinations of 
i.i.d. zero-mean processes, which are associated to time 
series modeling, do satisfy almost surely such a weak 
assumption. Our analysis 

• does not make any difference between non-stationary 
and stationary time series, 

• does not need the often tedious and cumbersome 
trend and seasonality decomposition (our trendlines 
include the seasonalities, if they exist). 

1.3 A model-free setting 

It should be clear that 

• a concrete time series cannot be "well" approximated 
in general by a solution of a "parsimonious" Eq. (1), 
i.e., a linear difference equation of low order; 

• the use of large order linear difference equations, or 
of nonlinear ones, might lead to a formidable compu- 
tational burden for their identifications without any 
clear-cut forecasting benefit. 

We adopt therefore the quite promising viewpoint of 
(Fliess & Join, 2008) where the control of "complex" sys- 
tems is achieved without trying a global identification but 
thanks to elementary models which are only valid during 
a short time interval and are continuously updated. We 
utilize here low order difference equations. Then the 
window size for the moving average (6) does not need to 
be "large". 

1.4 Content 

Sect. 2, which considers the identifiability of unknown 
parameters, extends to the discrete-time case a result in 
(Fliess, 2008). The convincing computer simulations in 
Sect. 3 are based on the exchange rates between US Dollars 
and €uros. Besides forecasting the trendhne, we predict 

• the position of the future rate w.r.t. the forecasted 
trendline, 

• the standard deviation w.r.t. the forecasted trendhne. 

Those results might lead to a new understanding of volatil- 
ity and risk management. ^ Sect. 4 concludes with a short 
discussion on the notion of trend. 

^ See the numerous examples and the references in (Fliess & Join, 
2008) for concrete illustrations. 

* Compare with (Markovsky, Willems, van Huffel, de Moor & 
Pintelon, 2005). 

^ See (Taleb, 1997) for a critical appraisal of the existing literature 
on this subject, which is of utmost importance in financial engineer- 
ing. (Extreme) risks are discussed in (Bouchaud & Potters, 1997; Da- 
corogna, Gengay, Miiller, Olsen & Pictet, 2001; Mandelbrot & Hud- 
son, 2004; Sornette, 2003) from quite different perspectives. It is the 
trendline which would exhibit abrupt changes in our setting (compare 
with the probabilistic standpoint; see, e.g., (Wilmott, 2006) and the 



2. PARAMETER IDENTIFICATION 

2. 1 Rational generating functions 

Consider again Eq. (1). The Z-transform X o{ x satisfies 
(see, e.g., (Doetsch, 1967; Jury, 1964)) 

- x'(0) - x{l)z-^ x{n- l)^-("-l^] , . 

a„_i2[X - x(0)] - a„X = ^ > 

It shows that X, which is called the generating function 
of X, is a rational function of z, i.e., X e R{z): 

P{z) 



X 



Q{z) 



(8) 



where 

P{z) = boz"-^ + + . . . + e R[z] 

Q{z) = z" cin-iz - a„ e M.[z\ 

Hence 

Proposition 2.1. x(t), t > 0, satisfies a linear difference 
equation (1) if, and only if, its generating function X is a 
rational function. 

It is obvious that the knowledge of P and Q permits to 
determine the initial conditions a;(0), . . . , x{n — 1). 

Remark 2.2. Consider the inhomogeneous linear difference 
equation 

x{t + n) — aix{t + n — 1) — • • • — a„x{t) 
= ^ m{t)a* + ^ m'{t) sm{ujt + (p) 

finite finite 

where m{t),w'{t) e R[t], a,uj,tp e R. Then the Z- 
transform X G M.{z) of x{t) is again rational. It is equiv- 
alent saying that x{t), t >Q, still satisfies a homogeneous 
difference equation. 

2.2 Parameter identifiability 



Generalities Let 



j? = Q(ai 



, an, bo, ... , bn~l) 



be the field generated over the field Q of rational num- 
bers by oi, . . . , a„, 60, . . . , bn-i, which are considered as 
unknown parameters and therefore in our algebraic setting 
as independent indeterminates (FHess & Sira- Ramirez, 
2003; Fliess, Join & Sira-Ramlrez, 2004; FHess & Sira- 
Ramlrez, 2008). Write ^ the algebraic closure of .^ (see, 
e.g., (Lang, 2002; Chambert-Loir, 2005)). Then X_ € i?(z), 
i.e., X is a rational function over A. Moreover A{z) is a 
differential field (see, e.g., (Chambert-Loir, 2005)) with 
respect to the derivation Its subfield of constants is 
the algebraically closed field R. 

Introduce the square Wronskian matrix M of order 2n + l 
(Chambert-Loir, 2005) where its ^''^-row, < x < 2?!, is 

(9) 



^z"X, ^X, ^z' 
dz^ dz^ dz^ 



dz^' 



-1 



It follows from Eq. (7) that the rank of ^A is 2n if, and 
only if, X does not satisfy a linear difference relation of 
order strictly less than n. Hence 



references therein). Our estimation techniques permit an efficient 
change-point detection (Mboup, Join & Fliess, 2008), which needs 
to be extended, if possible, to some kind of forecasting. 



Theorem 2.3. If x does not satisfy a linear difference 
equation of order strictly less than n, then the parameters 

ai, . . . , a„, bo, ... , bn-i 

are linearly identifiable. ^ 

Identifiability of the dynamics For identifying the dy- 
namics, i.e., oi, . . . , Gn, without having to determine the 
initial conditions consider the (n+l) x (n+ 1) Wronskian 
matrix M, where its ^*''-row, < /i < 71 + 1, is 

— 1 — Z X, .... \ — X 

It is obtained by taking the X-dependent entries in the 
n + l last rows of type (9), i.e., in disregarding the entries 
depending on 60, ... , bn-i. The rank of A/" is again n. Hence 
Corollary 2.4- oi, . . . , a„ are linearly identifiable. 

Identifiability of the numerator Assume now that the 
dynamics is known but not the numerator P in Eq. (8). 
We obtain bg, . . . , 6„_i from the first n rows (9). Hence 

Corollary 2.5. bo, ... , 6„_i are linearly identifiable. 

2.3 Some hints on the computer implementation 

We proceed as in (Fliess & Sira-Ramlrez, 2003; Fliess & 
Sira-Ramirez, 2008) and in (Fliess, Fuchshumer, Schoberl, 
Schlacher & Sira-Ramirez, 2008). The unknown linearly 
identifiable parameters are solutions of a matrix linear 
equation, the coefficients of which depend on x. Let 
us emphasize that we substitute to x its filtered value 
thanks to a discrete-time version of (Mboup, Join & 
Fliess, 2009).'^ 

3. EXAMPLE: FORECASTING 5 DAYS AHEAD THE 
$ - € EXCHANGE RATES 

We are utilizing data from the European Central Bank, 
depicted by the blue lines in the Figures 1 and 2, which 
summarize the 2400 last daily exchange rates between the 
US Dollars and the €uros. ® 

3.1 Forecasting the trendline 

In order to forecast the exchange rate 5 days ahead we 
apply the rules sketched in Sect. 2.3 and we utilize a linear 
difference equation (1) of order 3 (the filtered values of the 
exchange rates are given by the black lines in the Figures 1, 
2). Fig. 3 provides the estimated values of the coefficients 
of the difference equation. The results on the forecasted 
values of the exchange rates are depicted by the red lines in 
the Figures 1 and 2, which should be viewed as a predicted 
trendline. 



® It means following the terminology of (Fliess & Sira-Ramirez, 
2003; Fliess & Sira-Ramirez, 2008) that ai, . . . , a„, feoi ■ ■ • i ''n-l 
are uniquely determined by a system of 2n linear equations, the 
coefficients of which depend on -j^^ and ^^2™, < m < n — 1. 

See, e.g., (Gengay, Selguk & Hafner, 2002) for an excellent presen- 
tation of various filtering techniques in economics and finance. 
* The authors are perfectly aware that only computations dealing 
with high frequency data might be of practical value. This type of 
results will be presented elsewhere. 
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Figure 1. Exchange rates (blue -), filtered signal (black - 
-), forecasted signal (5 days ahead) (red -) 
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Figure 3. Parameter estimations ai (red —.—), a2 (blue 
) and 03 (black ) (5 days ahead) 

3.2 Above or under the predicted trendline? 

Consider again the "error" iy{t) in Eq. (2) and its moving 
average MA„,N{t) in Eq. (6). Forecasting MAy_7v(i) as in 
Sect. 3.1 tells us an expected position with respect to the 
forecasted trendline. The blue line of Fig. 4 displays the 
result for the window size N = 100. The meaning of the 
indicators A and V is clear. 

Table 1 compares for various window sizes the signs of the 
predicted values of MA^^N{t), which tells us if one should 
expect to be above or under the trendline, with the true 
positions of x{t) with respect to the trendline. The results 
are expressed via percentages. 

3.3 Predicted volatility w.r.t. the trendline 

Introduce the moving standard deviation 
MSTD^,]v(i) = 

V N + l 




Figure 2. Zoom of Figure 1 
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Figure 4. Predicted position w.r.t. the trendline (5 days 
ahead) - V: above, A: under 



Window's size 


Percentage 


50 


65.6% 


100 


88.3% 


200 


62.3% 


300 


67.1% 



Table 1. Comparison between the sign of the 
predicted value of MA^,Ar(t) and the true posi- 
tion of x{t) w.r.t. the trendline (5 days ahead). 

and forecast it as in Sect. 3.2. The results, which are 
displayed for a window size N = 100 in Table 2 and Fig. 
5 via the familiar confidence intervals, ^ confirm the time- 
dependence of the variance, i.e., the heteroscedasticity. 



Confidence intervals 


Prediction 


Real 


mean-3 x std,mean+3 x std 


99% 


98.7% 


mean-2 X std,mean+2 X std 


95% 


92.2% 


mean-std,mean+std 


68% 


64.4% 



Table 2. Confidence interval validations (5 days 
ahead) 
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Figure 5. Confidence interval (95%) (5 days ahead) 
3.4 Forcasting 10 days ahead 

Figures 6, 7, 8, 9 display the same type of results as 
in Sections 3.1, 3.2, 3.3 via similar computations for a 
forecasting 10 days ahead. The quality of the computer 
simulations only slightly deteriorates. 



^ There is of course no need for the underlying statistics to be 
Gaussian. Lack of space prevents us from exhibiting forecasts of 



quantities like skewness and kurtosis, which would be obtained by 
similar calculations. This will be done in some future publications. 
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Figure 7. Zoom of Figure 6 
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Figure 6. Exchange rates (blue -), filtered signal (black - 
-), forecasted signal (red -) (10 days ahead) 

4. CONCLUSION 

The existence of trends, which is 

• the key assumption in technical analysis, ^° 

• quite foreign, to the best of our knowledge, to the 
academic mathematical finance, where the paradigm 
of random walks is prevalent (see, e.g., (Wilmott, 
2006)), 



Trends in technical analysis should not be confused with what are 
called trends in the time series literature (see, e.g., (Gourieroux & 
Monfort, 1995; Hamilton, 1994)). 
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Figure 8. Predicted position w.r.t. the trendline (10 days 
ahead) - V: above. A: under 

is fundamental in our approach. A theoretical justification 
will appear soon (Fliess & Join, 2009). We hope it will 
lead to a sound foundation of technical analysis, which 
will bring as a byproduct easily implementable real-time 
computer programs. 



^^The existence of trends does not necessarily contradict a random 
character (see (Fliess & Join, 2009) for details). 
^■^See also (Dacorogna, Gengay, Miiller, Olsen & Pictet, 2001) for a 
most exciting study which employs high frequency data. There are 
also other types of attempts to put technical analysis on a firm basis 
(see, e.g., (Lo, Mamaysky & Wang, 2000)). See (Blanchet-Scalliet, 
Diop, Gibson, Talay & Tanre, 2007) for a comparison between 
technical analysis and model-based approaches with parametric 
uncertainties. 

Our technics already lead to such computer programs in automatic 
control and in signal processing. 
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Figure 9. Confidence interval (95%) (10 days ahead) 
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